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The  existing  body  of  literature  on  particle  transport  by  turbulent  flow  has  concentrated  on  the 
behavior  of  spherical  particles.  However,  the  dispersed  phase  in  several  suspensions  of 
industrial  interest  consists  of  nonspherical  particles,  in  particular,  flexible  slender  bodies  or 
threads.  The  present  work  considers  the  problem  of  disorientation  of  initially  aligned  material 
lines  in  isotropic  turbulence,  an  idealization  which  serves  to  test  a  general  model  for  particle 
orientation  effects  in  the  transport  of  threads  in  turbulent  flow.  Results  obtained  from  the 
model  agree  wq  1  with  data. 
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The  problems  of  orientation,  deformation,  and  motion 
of  flexible  and  slender  particles  in  turbulent  flow  arise  in 
areas  such  as  drag  reduction1,2  by  threads  of  concentrated 
polymer  solutions  or  by  suspended  fibers  and  in  various 
manufacturing  operations  involving  the  transport  of  threads 
by  turbulent  airflow.  In  the  present  article,  we  do  not  at¬ 
tempt  to  solve  this  problem  in  the  sense  of  obtaining  a  math¬ 
ematical  solution  in  closed  form  under  particular  circum¬ 
stances;  rather,  we  wish  to  obtain  a  closed  system  of 
equations  which  can  be  used  to  compute  various  statistical 
quantities  in  arbitrary  flow  situations.  This  requires  model¬ 
ing  of  various  terms.  If  carefully  constructed,  such  models 
can  be  more  general  than  solutions  resulting  from  particular 
assumptions  regarding  the  turbulence,  such  as  that  the  strain 
field  is  white,  or  Gaussian,  or  two  dimensional,  or  weak,  all 
of  which  would  permit  solutions  of  the  classical  form,  but 
none  of  which  is  particularly  realistic  when  applied  to  turbu¬ 
lence.  The  type  of  modeling  that  we  propose  avoids  making 
assumptions  of  this  type. 

In  this  article,  we  do  not  develop  a  general  model;  rath¬ 
er,  we  seek  a  situation  that  can  be  used  to  calibrate  an  equa¬ 
tion  for  the  evolution  of  particle  orientation,  which  is  one 
part  of  a  general  two-fluid  model  for  the  transport  of  thread¬ 
like  suspended  particles.  We  present  here  the  part  dealing 
with  particle  orientation,  specialized  for  a  particular  set  of 
circumstances.  The  specialization  corresponds  to  the  limit  of 
negligible  inertia,  perfect  flexibility,  and  perfect  extensibility 
of  the  particles.  In  this  limit,  the  problem  of  particle  orienta¬ 
tion  in  turbulent  flow  reduces  to  a  study  of  the  orientation  of 
material  lines  in  turbulent  flow.  The  effect  of  turbulence  on 
material  lines  and  surfaces  has  been  theoretically 1,4  and  ex¬ 
perimentally56  considered  in  previous  investigations.  In  par¬ 
ticular,  the  angular  dispersion  of  fluid  lines  has  been  experi¬ 
mentally  studied6  by  introducing  lines  of  hydrogen  bubbles 
in  the  nearly  isotropic  turbulence  behind  a  grid.  We  propose 
a  simple  model  for  the  dispersion  which  is  shown  to  be  con¬ 
sistent  with  the  experimental  results. 

II.  MODEL 

Let  us  denote  a  differential  element  of  the  fluid  line  by  ds 
which  has  components  dx, .  The  orientation  of  ds  is  given  by 
the  unit  vector  A  whose  components  are 


yjdXj  dx} 

We  need  to  derive  an  equation  for  the  evolution  of  A,  with 
time.  Such  an  equation  can  be  obtained  by  taking  the  sub¬ 
stantial  derivative  of  Eq.  ( 1 ).  After  differentiating  Eq.  ( 1 ), 
expressing  the  time  derivatives  of  dxt  and  dx}  dx}  in  terms  of 
the  spatial  derivatives  of  the  fluid  velocity  u,t  and  simplify¬ 
ing,  we  obtain 

~dT=a>  (2) 

which  if  rewritten  as 

DA, 

AjUtJ  =  +  At  {AJAkujk )  (3) 

can  be  seen  to  be  an  expression  for  the  relative  motion 
between  adjacent  points  on  ds  in  terms  of  the  rotation  of  ds 
and  a  stretching  along  ds.  It  should  be  noted  that  when  the 
two-fluid  model  is  used  to  describe  particle  transport  in  a 
carrier  fluid,  Eq.  (2)  is  the  equation  for  the  evolution  of 
particle  orientation. 

Equation  (2)  is  somewhat  complicated  by  the  require¬ 
ment  that  A,  be  a  unit  vector,  which  is  responsible  for  the  last 
term  in  the  equation.  We  could,  of  course,  work  with  the 
unnormalized  vector,  which  would  make  modeling  some¬ 
what  easier.  In  the  full  problem,  however,  the  normalized 
vector  A,  appears  very  often  in  the  momentum  equations. 
Using  the  unnormalized  vector  would  simplify  Eq.  (2)  and 
greatly  complicate  all  the  other  equations.  We  feel  that  the 
use  of  the  normalized  vector  produces  the  greatest  overall 
simplicity. 

Now  A,  and  u,  are  decomposed  into  their  means  A,  and 
and  U,  and  their  fluctuations  A  '  and  u\  as  follows: 

A.  =  A,  +  /t (4) 

ut  =  U,  +  u'r  (5) 

On  squaring  and  averaging  Eq.  (4),  we  obtain 

1  =  A,A,  +  AjAj.  (6) 

Thus,  although  A,  is  a  unit  vector,  Af  is  not  a  unit  vector. 

In  the  case  of  the  problem  when  both  the  turbulence  and 
the  suspended  matter  are  homogeneous  and  there  is  no  mean 
shear,  the  averaged  form  of  Eq.  (2)  reduces  to 
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d,(A, )  =  [  A  ;u'tJ  -  u]A  '  -  A,  A,  AluJ'i 

-  a,  A*  x;w;A  -  a,  a*  ] 

-  (  itty'*  )  -  {A,  A'jA'ku'/k 

-+■  Ay  A  jA  k  Uj  k  -h  A*  A  tA  j  ujk  }.  ( 7 ) 

The  third-order  correlations  in  curly  braces  in  Eq.  (7)  may 
be  dropped  if  we  assume  that  the  fluid  line  statistics  are  joint 
normal.  We  assume  the  carrier  fluid  turbulence  to  be  iso¬ 
tropic;  hence  we  expect  the  remaining  correlations  in  Eq. 

( 7 )  to  be  axisymmetric  around  A, .  We  also  assume  that  the 
fluid  lines  are  initially  aligned  in  a  particular  direction;  hence 
all  the  correlations  must  be  0  at  t  —  0.  Eventually  the  turbu¬ 
lence  completely  destroys  the  orientation  of  the  fluid  lines 
and  therefore  the  correlations  must  also  tend  to  0  as  /-*  oc . 
From  Eq.  (6)  it  can  be  seen  that  the  magnitude  of  A  '  is 
( 1  —  A: ) 1  /2,  where  A2  denotes  A,  A, .  We  normalize  by  the 
appropriate  power  of  the  magnitude  so  that  the  correlations 
will  behave  correctly  when  A  =  1,  which  occurs  at  /  =  0, 
and  when  A  — 0,  which  occurs  as/-*  oc.  Furthermore,  in  Eq. 
(7)  we  assume  that  there  is  a  time  scale  associated  with  the 
second-order  correlations  in  square  brackets  and  another 
time  scale  corresponding  to  the  third-order  correlation  in 
parentheses. 

The  above  considerations  lead  to  the  following  models 
for  the  terms  on  the  right-hand  side  of  Eq.  (7): 

[•]  =  -  A,  ( 1  -A2)1'2/,  (8) 

where  A:  =  Ay  Ay  and  f  =  inverse  time  scale  of  [  •  ] ; 

(  )  =  -  A,  ( 1  -  A2)3/V  (9) 

where  n  is  the  ratio  of  the  inverse  time  scales  of  (  )  and  ( •  ] ; 
{•}  =  0.  (10) 

After  substituting  Eqs.  (8 >-(10)  into  Eq.  (7)  and  forming 
an  evolution  equation  for  A2,  we  obtain 

d,{A2)  =  —  2A2(  1  -  A2)l/2/[l -n(\ - A2)].  (11) 

In  Eq.  (11),  since  d, (A2)=0  at  /  =  0,  we  require 
d2(A2)  <0  at  /  =  0  in  order  that  A2  decay  with  time.  This 
constraint  is  assured  by  a  value  of  {  for  the  fractional 
exponent  in  Eq.  ( 1 1 ) ;  no  other  value  would  satisfy  this  con¬ 
straint. 

Since  in  our  simple  problem  we  expect /  to  scale  with  the 
inverse  time  scale  of  the  fluid  turbulence,  we  take 

f-ae/q %  (12) 

where  £  is  the  isotropic  dissipation  rate  and  q2  is  twice  the 
turbulent  kinetic  energy.  The  expression  for  the  inertial  fre¬ 
quency  scale  of  the  decaying  turbulence  behind  the  grid  in 
the  modeled  experiment6  is 

e/q2  =  0.625/(2  +  /),  (13) 

where  /  is  the  time  elapsed  in  seconds  after  the  introduction 
of  the  material  lines  into  the  flow. 

Substituting  Eqs.  (12)  and  (13)  into  Eq.  ( 1 1 ) ,  we  ob¬ 
tain 

-  125«A2(1  —  A2)l/J[  1  — /»(!  —  A2)] 


(14) 


Equation  ( 14),  along  with  the  adjustable  parameters  a 
and  n%  is  an  equation  that  describes  the  progressive  disorien¬ 
tation  of  the  fluid  lines  with  time.  The  model  parameters  a 
and  n  are  chosen  so  that  the  behavior  of  A2(/)  obtained  by 
solving  Eq.  ( 14)  matches  with  the  experimental  data.6  The 
initial  condition  for  Eq.  (14)  is 

A2(/  =  0)  =  1.  (15) 

III.  RESULTS 

In  the  experimental  setup,6  the  flow  was  in  the  one-di¬ 
rection,  with  the  lines  introduced  initially  in  the  two-direc- 
tion,  as  shown  in  Fig.  1 .  The  results  of  that  experiment  were 
presented  as  a  plot  of  (0 ,  )2  (/  *),  where  {8^  )2  is  the  vari¬ 
ance  of  8y  (see  Fig.  1 ).  The  independent  variable  in  the  ex¬ 
perimental  plot  was  /  *,  which  represents  time  normalized  in 
a  manner  which  accounts  for  the  decay  of  the  grid  turbu¬ 
lence  and  is  given  by 

/*=  181  [0.707  —  ( 16  -f  8/ )  0 125 ] .  (16) 

Equation  ( 14),  along  with  its  initial  condition  ( 15),  was 
solved  analytically  for  A2(/).  In  order  to  compare  the  model 
with  experimental  data  we  need  an  expression  for  (8'}  )2  in 
terms  of  A2.  We  can  write 

#  ;=  —  arcsin  (A  \/yj\  —  A  J2 ) ,  (17) 

where  A  J  and  A  ;  are  the  fluctuations  in  the  components  of 
the  fluid  orientation  vector  along  the  one  and  three  axes, 
respectively. 

After  squaring  Eq.  (17);  expanding  its  right-hand  side 
in  a  Taylor  series  up  to  fourth  order  in  A  [  and  A  { ,  assuming 
that  the  distributions  of  A  ;  and  A  3  are  statistically  indepen¬ 
dent;  and  averaging,  we  obtain 

7^17-  Tp  +  IT2  Ip  +  Jf/3.  (18) 

In  order  to  express  ( 0 , ) 5  in  terms  of  A2  we  now  need  ex- 
pressions  for  the  correlations  on  the  right-hand  side  of  Eq. 
( 1 8 )  in  terms  of  the  mean  orientation  vector  A, .  An  equation 
for  A  ■  A  '  can  be  derived  as  follows.  In  the  homogenous 
problem  which  is  considered  in  the  present  work,  the  anisot¬ 
ropy  in  A  '  A  j  is  due  to  the  anisotropy  in  A,  A, ;  hence 

JJJJ  =(1  -A2)[So/3+0(AtAj  -A2<$,/3)]«  09) 

To  compute  /?,  consider  almost  aligned  fibers;  A, 
=  ( 1  —  e,0,0),  where  £  is  a  small  parameter.  In  such  a  case 


FIG.  1,  Notation  for  the  fluid  line  orientation. 
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FIG  2.  The  transient  behavior  of  angular  dispersion. 


Eq.  (19)  gives  a  diagonal  form  for  A  '  Aj,  which  to  first 
order  in  e  is 


A’,2  =  2  +  $/?),  A  22  =  2e(J  —  }/?), 

Aj2  =  2e{\-\P). 


(20) 


We  require  that  /?<  1  in  order  that  A  j2  and  A  '2  be  always 
non-negative.  It  seems  plausible  that  /?  is  a  weak  function  of 
A hence  we  assume  that  0  is  a  constant  and  set  it  equal  to  1 . 
Taking/?  =  1  in  Eq.  ( 19)  results  in 

A ;  a;  =  ( i  -  a2>[V3  +  (a, a,  <21) 

In  the  present  case,  due  to  symmetry,  A ,  =  0  and  A3  =  0  for 
ail  time.  Therefore,  from  Eq.  v  2 1 ), 


A\2=  A?  =  (  1  -  A2)2/3. 


(22) 


If  we  assume  that  A  J  has  a  Gaussian  distribution,  then 


a;4  =  3(  I7V.  (23) 

Substituting  Eqs.  (22)  and  (23)  into  Eq.  ( 18),  we  obtain 

(  -  A:)73  +  2(l  -  A:)4/9.  (24) 


Numerical  results  of  {6  3  )2{t  *)  were  obtained  for  var¬ 
ious  values  of  the  parameters  a  and  n.  The  choice  of 

a  ~  9.4  and  n  =  0.5  (25) 

gives  a  good  fit  between  the  prediction  of  the  model  and  the 
experimental  data,  as  shown  in  Fig.  2.  It  should  be  noted  that 
Eq.  (24)  is  obtained  by  a  Taylor  expansion  of  Eq.  (17)  for 
small  A  ;  and  therefore  we  do  not  expect  a  good  match  for 
large  time.  There  is  no  reason,  however,  to  expect  Eq.  (11) 
to  be  any  less  accurate  for  large  time.  Substituting  Eqs.  (12) 
and  (25)  into  Eq.  ( 1 1 ),  we  obtain 

<?,(A2)  =  -  I8.8e A2(l  -  A2) ,/2[  1  -0.5(1  -  A2)]/?2. 

(26) 

In  conclusion,  we  have  been  able  to  derive  an  ordinary 
differential  equation,  Eq.  ( 26 ) ,  for  the  mean  orientation  vec¬ 
tor  A,,  which  characterizes  the  disorienting  effect  of  turbu¬ 
lence  on  initially  aligned  material  lines  in  an  isotropic  turbu¬ 
lent  fluid.  Perhaps  more  valuable,  we  have  used  this  special 
case  to  calibrate  a  model  which  may  have  more  general 
utility. 
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